Comment on "Dielectric behavior of paraelectric KTaOs, CaTiOa, and 
(Lni/2Nai/2)Ti03 under a dc electric field" 
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Chen Ang, Bhalla and Cross [Phys. Rev. B 64, 184104 (2001)] have studied the low-frequency 
(20 Hz - 100 kHz) dielectric dispersion of the KTaOs crystal under a dc electric field. Perform- 
ing fits of the electric field dependence they came to conclusion that an appreciable contribution 
to the dielectric permittivity originates from polar clusters. In this Comment we show that the 
dielectric permittivity at low frequencies (100 Hz - 1 MHz) equals to that in the THz region, 
close below the polar phonon response. This excludes the possibility of any appreciable dielectric 
dispersion due to polar clusters. In addition, we demonstrate that correct treatment using Landau- 
Ginzburg-Devonshire theory allows to fit the electric field dependence of the dielectric constant 
without assuming any polarization mechanism besides the polar phonon modes. 

PACS numbers: 77.84.Dy, 77.80.Bh 



The family of incipient ferroelectric materials with per- 
ovskite structure has attracted a great attention due 
to its high value of dielectric constant, low dielectric 
losses and high degree of tunability by the external elec- 
tric field and temperature. The best known members 
of this family are SrTiOg (STO), KTaOg (KTO)i and 
CaTiOa (CTO)^ which are also called "quantum para- 
electrics". Microwave measurements of these materials 
by Rupprecht and Bell2i demonstrated that there is no 
sign of appreciable dielectric dispersion up to the GHz 
range. Therefore it is considered that dielectric proper- 
ties of pure incipient ferroelectrics are determined by the 
contribution of the lowest polar-phonon mode (the soft- 
mode) whose frequency is decreasing on cooling roughly 
according to the Cochran law. Recently, combining the 
far infrared (FIR) spectroscopy with low-frequency (LF) 
dielectric measurements some of present authors con- 
firmed the absence of appreciable dielectric dispersion be- 
low the soft-mode response in STO ceramics** and CTO 
single crystals^. 

Chen Ang, Bhalla and Cross in their paper— in- 
vestigated the dielectric response of KTO, CTO and 
(Lni/2Nai/2)Ti03 (Ln = La and Nd) (LNTO) under a 
strong dc electric field (up to 60 kV/cm). To describe the 
field dependence of the permittivity, they attempted to 
use an equation which is a simple polynomial expansion 
by even powers of electric field up to the fourth order. 
They were able to describe the field dependence in CTO 
and LNTO, but failed to perform a satisfactory fit for 
KTO. This fact brought the authors to conclusion that 
there should be another Langevin-type term describing 
the relaxational contribution of polar clusters to the total 
polarization. The final equation was of the form^ 



e{E) = ei - e2E'^ + e^E^ + [PrX / e{Q)][coiih{Ex)]-'^ (1) 

where x — PrL^E/2kBT, Pr is the effective polariza- 
tion of the polar clusters with the size L and ks is the 
Boltzmann's constant. Using Eq. n]to fit the measured 



field dependence of the permittivity, they obtained the 
following parameters of polar clusters: Pr = 0.4 /iC/cm^, 
L = 4.7 - 6 nm and Aeduster = 656 (at 14 K). This would 
mean that 23% of the total dielectric response at low fre- 
quencies are due to some relaxational dispersion of polar 
clusters below the soft-mode response. 

In the following we will show that our recent FIR and 
LF dielectric measurements of KTO exclude the possibil- 
ity of any appreciable contribution to the LF permittivity 
besides the contribution of the soft phonon mode. More- 
over we will demonstrate that a correct treatment using 
the Landau-Ginzburg-Devonshire thermodynamic theory 
allows to describe the field dependence of the dielectric 
constant in KTO without using additional Langevin-type 
terms. 

Our measurements have been performed on nominally 
pure KTO single crystals^. The samples were prepared 
in the form of thin < 001 >-oriented plane-parallel plates 
with dimensions of 7 x 7 x 0.2 mm"^. LF dielectric mea- 
surements were carried out using a HP 4192A impedance 
analyzer in the frequency range 100 Hz - 1 MHz with gold 
electrodes deposited on the main surfaces of the sam- 
ple. The FIR refiectivity was measured with Bruker IFS 
113v Fourier Transform spectrometer in the range 30 - 
650 cm~^. The complex dielectric function in the fre- 
quency range 5-30 cm~* was obtained by means of 
time-domain THz transmission spectroscopy on a 65 pm 
thick plane-parallel sample. The frequencies of the three 
polar phonon modes and their contribution to the di- 
electric function were obtained from the fit to the FIR 
reflectivity using the factorized form of the model dielec- 
tric function combined with the directly measured values 
of the THz complex permittivitj^ . The temperature de- 
pendence of the soft-mode frequency was found to be in 
perfect agreement with the hyper-Raman data—. 

In Fig.nthe measured permittivity at 43.5 kHz is com- 
pared with the low-frequency limit of our fitted FIR di- 
electric function. The permittivity data from Ref. [Q are 
also added for comparison. One can observe good agree- 
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FIG. 1: Permittivity of KTaOs crystal at 43.5 kHz (full line), 
LF permittivity data taken from Ref. .fi (dotted line) and con- 
tribution of the polar-phonon modes (open circles). 
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FIG. 2: Bias electric-field dependence of the KTaOs permit- 
tivity at 14 K taken from Ref. (the first run). Full line 
corresponds to the fit using Eq.Eland0] 



ment between LF permittivities of KTO crystals and be- 
tween LF and FIR permittivities in the whole measured 
temperature range. Somewhat higher value of FIR per- 
mittivity compared to the LF values at 20 K lies within 
the instrumental errors. Anyway, Fig. ^ clearly shows 
that at all measured temperatures there is no space for 
any appreciable dispersion below the soft-mode response 
down to kHz range. 

To demonstrate the fitting of the field dependence of 
permittivity, following Fleury and Worlock^*^ let us ex- 
pand the crystal free energy per unit volume in a power 
series of lattice polarization: 

G = Go + lx{T){P^ + P'y + P!) + l^P' + Py+ Pz) 

+ yiP'Py + PyP! + P'P!) + lap! + Pi + Pt) 

+ other sixth-order terms + . . .(2) 

In our case polarization has only one component di- 
rected along the z-axis (Pa; = Pj, = 0). Thus we can 
obtain relation between the electric field and polariza- 
tion: 

E,, = {dG/dP,)T = X{T)P. +^P! + CP! + ... (3) 

which is a special case of the general relation for arbi- 
trary polarization^"' Then one can derive the expres- 
sion for the inverse susceptibility: 

Lht!e)-1] )^^ - ^^^/^^^ = X(T)+3CP.^+5CP^+. . 

(4) 

Here x(^) = l/eo[e(T, 0) — 1] is the inverse susceptibil- 
ity in the zero bias field and eo = 8.85 • 10~^^ F/m is the 
permittivity of vacuum. To calculate the permittivity un- 
der dc bias field, one has to use Eq.^to calculate the po- 
larization and then insert it into Eq. ^ The calculations 



should be performed numerically because Eq.|3|is a non- 
linear polynomial equation of the 5-th order. Therefore 
the analytical expression given in Ref. cannot be valid 
for highly anharmonic crystals and can be applied only 
in the weak field approximation. This approximation ap- 
parently holds in the case of CTO and LNTO because of 
their weak nonlinearity for bias fields i? < 60 kV/cm. 
However, in the case of KTO the lattice anharmonicity 
is essential. The fit using Eq. Oand^lof KTO permittiv- 
ity under the dc electric field at 14 K taken from Ref. 
(the first run of measurements) is shown in Fig.|21 The fit 
uses two independent parameters ^ — 10^° and C = 10^"'^ 
in SI units and the value e(14 K, 0) — 2821 taken from 
Ref. 0. The fitted nonlinear coefficients are in reason- 
able ^reement with the previously published data (see 
Ref. 10 and references therein). Fig.[21demonstrates that 
the expansion up to the 6-th order in Eq. |21 is sufficient 
to ensure a good agreement between the theory and ex- 
periment. 

The critical limit for the weak field approximation can 
be found from the condition that 4-th and 6-th order 
terms in Eq.[21are small compared to the quadratic term 

^P^ + CPi = Ax, (5) 

where A ^ 1 is a dimensionless constant. For the 
typical values of nonlinear coefficients^^ we can neglect 
the 4-th order term in Eq. and using Eq. |2| we obtain 




Using our fit parameters for KTO and A = 0.1 we 
calculate Ecr — 1.3 kV/cm. One can see that this value 
is much smaller than the highest dc electric field in Ref. 
(60 kV/cm). 

It is believed that polar nano-clusters are induced by 
unavoidable dipolar impurities in nominally pure KTOi^. 
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The presence of such clusters appears in the LF dielectric 
spectra of KTO as a universal low-temperature dielec- 
tric dispersion^ with thermally activated behavior. Be- 
low 30 K this dielectric dispersion shifts below the fre- 
quency range of conventional LF dielectric spectroscopy. 
Dielectric spectroscopy can only establish that the con- 
tribution to permittivity of polar clusters is very small 
comparing to the phonon contribution above 30 K and 
cannot be observed in LF dielectric spectra below 30 K. 
On the other hand, the polar clusters can be detected 
by light scattering as it has been done in Raman^'^ and 
hyper-Raman^^ scattering experiments. These experi- 
ments, however, bring no information about the dynam- 
ics of the clusters because their relaxation is quasi-static 
comparing to the frequencies resolved in the light scat- 
tering experiments. The local lattice symmetry breaking 
due to the polar clusters relaxes the selection rules and 
induces activation of new lines. These changes are in 
agreement with existence of the polar clusters and do 
not contradict with our dielectric spectroscopy results. 

Taking into account all the mentioned arguments, we 
can conclude that the estimated size of polar clusters 



found in Ref. ^ does not have any reliable physical back- 
ground and its possible correspondence with the Raman 
resultsi^ is accidental. 

Analysis of the field dependence of permittivity us- 
ing equation similar to Eq. ^ was previously applied 
also to Sri^xCaxTiOa system.^"'- where, due to the Ca- 
doping, appreciable dielectric dispersion indeed exists. 
However, the applied electric fields were much smaller 
(E < 1 kV/cm) than those in Ref. Estimation 
of the critical electric field from Eq. El gives for STO 
Ecr — 0.2 kV/cm (nonlinear parameters are taken form 
Ref. 113). Thus even 1 kV/cm field can introduce some 
error if one uses Eq. ^ In our opinion, the correct ther- 
modynamic approach described here is always preferable 
for the analysis of nonlinear behavior under strong elec- 
tric fields. 
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